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SUMMARY 


In order to obtain a broader understanding of the factors which 
affect spin and recovery motions, a technique has been devised for deter- 
time histories of the attitudes and velocities of free— spinning- 
tunnel models from film taken by two motion-picture cameras operating 
simultaneously. This paper presents the method devised and the results 
of an initial application of the method in which time histories of atti- 
tudes and velocities have been determined for one medium-attitude, mod- 
erately oscillatory, developed spin and the recovery therefrom for a 
model representative of a contemporary fighter airplane. 

The time-history curves indicated that the oscillatory motion of the 
spin was not completely regular, inasmuch as the period and amplitude of 
any one cycle of the time histories were not exactly the same as those for 
the preceding and following cycles. After rudder reversal from with to 
against the spin, some increases occurred in the amplitude of the oscilla- 
tions, particularly in the lateral attitude angles of sideslip and wing 

tilt. The recovery was completed about 1^ turns after the rudder was 

reversed. One factor which appeared to have an important part in the 
attainment of the recovery from the spin in this investigation was a 
relatively high inward sideslip reached during the oscillations following 
rudder reversal. Analysis indicates that the high inward sideslip appar- 
ently caused a relatively large negative aerodynamic rolling moment which 
caused a large decrease in rolling velocity; as a result, there was a 
corresponding large decrease in the gyrodynamic pitching moment which had 
been holding the model in a nose-up spinning attitude. 
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INTRODUCTION 


Many investigations have heen made in the Langley 20-foot free- 
spinning tunnel in which airplane spin and recovery characteristics have 
heen determined experimentally with dynamically scaled models. The results 
of a few such investigations are presented in references 1 to 5* From the 
results of the free-spinning-tunnel tests , it has heen possible to evaluate 
some of the factors affecting spins and recoveries. However, with the 
advent of modern high-speed airplanes with design features such as extreme 
fuselage -heavy loadings, long nose lengths, and thin swepthack wings, the 
complexity of the spinning problem has heen accentuated. As one means of 
obtaining a broader understanding of the factors affecting spins and recov- 
eries, a technique has been devised for determining time histories of the 
attitudes and velocities of the models from film taken by two simultane- 
ously operating mot ion -picture cameras. As an initial application of this 
technique, an investigation has been made to determine time histories of 
the attitudes and velocities during one spin and the recovery therefrom 
for a model representative of a contemporary fighter airplane. 

In this paper are presented descriptions of the camera setups and 
associated equipment, the equations and methods used in determining the 
model attitudes and motions from the film, and time-history curves of the 
variables of the spin and recovery investigated. The attitudes and veloc- 
ities are plotted with respect to airplane body axes, the earth, the rel- 
ative wind, and the axis of the resultant rotary motion. The curves show 
the attitudes and velocities of the model for a continuous time period of 
about 19 seconds (full-scale time), of which the first 12.5 seconds or 

so represent about 5^ turns of a developed spin and the remaining time 

represents a 1—turn recovery from the spin. Discussions of the nature 

of the spin and recovery obtained and of factors ostensibly affecting 
these motions are included. 


SYMBOLS 


The time-history plots are presented with respect to 
body system of axes and to space and wind attitude angles 


the airplane 
(fig. 1). 


X,Y,Z 


c 


l 



airplane body axes 
rolling -moment coefficient, 


NACA TN 3611 


3 


*-111 


M 


pitching -moment coefficient, — 

|pV R 2 Sb 




yawing -moment coefficient, 


N 


ioVj, 2 ®. 


longitudinal-force coefficient. 


|pV r 2 s 


°Y 


lateral-force coefficient. 


|pV R 2 S 


r n 

normal-force coefficient, — 

|pv R 2 s 


L rolling moment, ft -lb 

M pitching moment, ft-lb 

N yawing moment, ft-lb 

Fy longitudinal force, lb 

Fy lateral force, lb 

F z normal force, lb 

S wing area, sq ft 

b wing span, ft 

p air density, slugs/cu ft 

u,v,w components of resultant linear velocity V R along the X, Y, 

and Z body axes, respectively, positive in directions of 
positive X-, Y-, and Z-forces, ft/sec 

V R resultant linear velocity, ft/sec 
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components of resultant angular velocity ft about X, Y, 
and Z body axes, respectively, positive in directions of 
positive L-, M-, and N-moments, radians/ sec 

a 

resultant angular velocity (if axis of resultant rotation 
is vertical, ft = ^ e ), radians/sec 


airplane relative -density coefficient, 

m 

mass of airplane, Weighty s ]_ U g S 

0 

c 

mean aerodynamic chord, ft 

x/c 

ratio of distance of center of gravity rearward of leading 
edge of mean aerodynamic chord to mean aerodynamic chord 

z/c 

ratio of distance between center of gravity and X body axis 
to mean aerodynamic chord (positive when center of gravity 
is below X body axis) 

I X ,I Y' I Z 

moments of inertia about X, Y, and Z body axes, respec- 
tively, slug-ft^ 

IS) 

2 

product of inertia about X and Z body axes, slug-ft 

I X " I Y 
mb^ 

I Y - I Z 
mb^ 

inertia yawing-moment parameter 
inertia rolling-moment parameter 

h-h 

mb^ 

inertia pitching-moment parameter 

g 

acceleration due to gravity, taken as 32.17 ft/sec^ 

6 e 

projection in vertical plane of total angular deflection 
of X body axis from a reference position in horizontal 
plane, positive when sense of deflection is airplane nose 
upward from reference position, deg or radians as indi- 
cated (fig. 1(a)) 
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t e 


0e 


0 

“1 

a l' 


a 


P 


T 


VR 




projection in horizontal plane of total angular defelction 
of X body axis from reference position in horizontal 
plane, positive when clockwise as viewed vertically from 
above airplane, deg or radians as indicated (fig. 1(b)) 

angular deflection in YZ body plane of Y body axis from 
horizontal plane, positive when clockwise as viewed from 
rear of airplane (if X body axis is vertical, 0 e is 

measured in horizontal plane from a reference line in this 
plane), deg or radians as indicated (fig. 1(c)) 

angle between Y body axis and horizontal, positive for erect 
spins when airplane right wing is down and for inverted 
spins when airplane left wing is down, deg 

acute angle between X body axis and vertical, deg 
acute angle between Z body axis and horizontal, deg 


angle between X body axis and axis of resultant rotation, 
,-l P 


cos -*• or sin 


. -1 jx 2 + q 2 

m - 1 - -» — 


n 


; deg 


angle of attack between relative wind Vp projected into 

XZ plane of symmetry and X body axis, positive when 
relative wind comes from below XY body plane, deg 

angle of sideslip at center of gravity (angle between rela- 
tive wind V R and XZ plane of symmetry), positive when 

relative wind comes from right side of plane of symme- 
try, deg 

angle of inclination of relative wind V R from XY body 
plane, positive when inclined below XY body plane, 

sin -1 deg 
V R 


angle between relative wind projected into XY body 

±\ 

plane and X body axis, positive when comes from 

right side of plane of symmetry, sin - - 1 - v — , deg 

+ v 2 


MCA TO 5611 


angle of inclination of axis of resultant rotation from 
XY body plane , positive when inclined below XY body 

plane, sin -1 £, deg 


angle between axis of resultant rotation projected into 
XY body plane and X body axis, positive when vector of 
axis of resultant rotation lies on right side of plane of 


symmetry, 


sin 


i/p 2 + q 2 


deg 


time, sec 

acute angle between direction of resultant linear velocity 
and axis of resultant rotation, deg 

velocity of airplane center of gravity along axis of result- 
ant rotation, V R cos cr, ft/sec 

angle between axis of resultant rotation and vertical, deg 
angle between resultant linear velocity and vertical, deg 


vertical component of velocity of airplane center of 
gravity , f t / s e c 

horizontal reference plane at eye level in tunnel, 23.48 feet 
above center of bottom-camera lens (fig. 19) 

fixed reference point where a vertical line through the side 
camera stanchion and pivot point intersects Q is 

10.95 feet from center of tunnel (fig. 19) 

vertical reference plane bisecting tunnel and bottom-camera 
lens, parallel to plane of film of side camera when side 
camera is undeflected (that is, when side camera is aimed 
at center line of tunnel in plane Pjj, see fig. 19) 

horizontal projection of distance from Q to side-camera 
lens, positive when lens is toward tunnel (fig. 17), ft 

vertical projection of distance between line of sight of 

side camera to point on model or wall and Pg, positive when 
lens is deflected above p H (fig. 17), ft 
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Tj deflection from horizontal of line from side-camera lens to 

a specific point on model, positive when above % 

(figs. l6 and 19), deg 

A deflection from Py of plane which intersects Py in a 

horizontal line passing through center of bottom- camera 
lens and which contains a specific point on model, posi- 
tive when deflected toward side camera (figs. l8 and 19), 
deg 

e horizontal projection of angle between a line from Q to 

a specific point on model and a vertical reference plane 
perpendicular to Py (figs. 18 and 19), deg 

Sy vertical distance between a specific point on model and P H' 

positive when point on model is above P H (fig. 19 )> ft 

s„ horizontal projection of distance between a specific point 

on model and Q (fig. 19 ), ft 

s' projection of s.^ in direction perpendicular to P Tr 

n n V 

(fig. 19), ft 

Pjj radius vector of polar coordinate system locating center of 

gravity of model in horizontal plane from center line of 
tunnel, ft 

Pjj as seen on film-magnifier screen (fig. 18), deg 

£ angle of polar coordinate system locating center of gravity 

of model in horizontal plane in tunnel (fig. 18), deg 

v acute angle between projection of X body axis in a 

horizontal plane and p^, positive when measured to right 

of projection of X body axis as seen from a position 
directly below nose end of X body axis (figs. 18 and 21), 
deg 
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B angle in horizontal plane through center of gravity, between 

intersection of vertical plane through X body axis with 
horizontal plane and intersection of plane of symmetry with 
horizontal plane; also, angle in horizontal plane through 
center of gravity, between intersection of vertical plane 
through Y body axis with horizontal plane and intersection 
of YZ body plane with horizontal plane, positive when 0 
is positive (figs. 21 and 22), deg 

A angle in horizontal plane through center of gravity, between 

intersection of vertical plane through Z body axis with 
horizontal plane and intersection of plane of symmetry with 
horizontal plane, positive when 0 is positive (figs. 21 
and 22 ) , deg 


Subscripts: 


n 

t 

wr 
w l 
top 
H 
eg 


direction of nose end of X body axis 
direction of tail end of X body axis 
direction of right-wing end of Y body axis 
direction of left-wing end of Y body axis 
direction of upper end of Z body axis 
projection in horizontal plane 
center of gravity 


A dot over a symbol represents a derivative with respect to time; 

for example, p = — . 

dt 

The attitude angles 0 e , \|r e , and 0 e ( Euler ean-type space angles) 

as defined and used herein are intended to provide a means of inter- 
preting any changes from the reference attitudes for either an initially 
erect or an initially inverted spin. Under this general system, the ref- 
erence attitude for 0 e is arbitrarily taken as zero for either an erect 

or an inverted spin; the model or airplane would go from an erect to an 
inverted attitude or vice versa if 0 e becomes greater negatively 

than - 90 ° or if 0 e exceeds ±90°. 
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APPARATUS AND METHODS 
Model 


A scale model which was representative of a contemporary fighter 

24 

design with a 35 ° sweptback wing was used in the investigation. A three- 
view drawing and a photograph of the model are shown in figures 2 and 3- 
Dimensional characteristics of the airplane represented are listed in 
table I. The model was made of balsa, spruce, and aluminum. Lead bal- 
last was added in such a manner that the model was considered to repre- 
sent dynamically a corresponding full-scale airplane at an altitude of 
15^000 feet. The mass characteristics of the model as ballasted are 
presented in table II in terms of full-scale values. A remote -control 
mechanism was installed in the model to actuate the model controls as 
desired. 


Wind Tunnel and Tests 

The model spin and recovery were performed in the Langley 20-foot 
free-spinning tunnel. This tunnel and its operation are generally sim- 
ilar to that described in reference 1 for the Langley 15-foot free- 
spinning tunnel, the only difference in operation being that the models 
are launched by hand with rotation instead of being launched mechanically 
into the vertically rising airstream. The models are launched with the 
controls set so as to maintain a spin and the tunnel airspeed is adjusted 
by an operator so that the weight of the model is supported. The model 
motion evolves into a filly developed spin or ceases rotating and dives 
in a "no spin." When a spin is obtained, it is observed for a number of 
turns, after which an attempt is made to effect a recovery from the spin 
by moving one or more of the model controls to oppose the spinning motion. 
Angles of attack for which spins are obtained range from about 70° for 
flat spins to about 30° for steep spins. 

Free-spinning-tunnel investigations have shown that, in general, 
developed spins may be classified as "steady" or "oscillatory" (ref. 2) 
depending on whether the model appears to be spinning smoothly or oscil- 
lating in roll, yaw, and pitch while spinning. The investigations have 
also shown that the motion following control movement for recovery may 
range from one in which the model, immediately after control movement, 
decreases its rate of rotation and noses down rapidly for a recovery to 
one in which movement of the controls has little or no effect on the 
spinning motion. In some instances, after controls are moved for recov- 
ery, oscillations may occur in the attitude of the model and in its rate 
of rotation with an overall trend toward a decrease in both rate of 
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rotation and angle of attack which eventually evolves into a steep low- 

angle -of -attack pull-out dive without rotation (\|r becomes zero). 

Recovery, however, is sometimes considered as having been achieved if the 
model appears to have entered a glide or a roll at an angle of attack 

below the stall, even though \| t e has not yet been reduced to zero. In 

some cases, recovery from the spin may be associated with a quick rolling 
or pitching motion to an inverted attitude; and recovery is considered to 
have been achieved from the original erect spin regardless of whether ljr 

has changed its direction in space. The motion of a specific airplane 
during a recovery is apparently dependent on the spin motion before the 
controls are moved, on the particular control manipulation employed for 
recovery, and on the effectiveness of these controls (refs. 3 to 5). 

For the present investigation, it was decided that a spin and a 
recovery of a type which has been very common in the past would be stud- 
ied. This was a spin of medium attitude and moderate oscillations, and 
the recovery was the type in which the model pitched down to a low angle 
of attack and ceased rotating within a couple of turns after the control 
was moved in attempting the recovery. Control settings used on the model 
for the developed spin were elevator neutral, ailerons neutral, and rud- 
der 30° with the spin (corresponding to right rudder pedal forward in a 
spin turning toward the right-hand side of the airplane pilot). Control 
manipulation used for recovery was reversal of the rudder to 15 ° against 
the spin. 


Camera Setup and Determination of Motion 

The technique by which the time histories of the attitudes and 
velocities of the model during the spin and recovery were determined is 
explained in detail in the appendix. In brief, the technique is as 
follows: Motion pictures were made of the model spin and recovery test 

with two cameras operating simultaneously. One camera was mounted at 
the bottom of the vertical tunnel and pointed upward along the center 
line of the tunnel; the other camera was at one side of the tunnel and 
about on a level with the spinning model. Time histories of the space 
coordinates of the model in the tunnel were obtained by means of timers 
mounted in view of each camera, suitable reference lines painted on the 
tunnel walls, and suitable attitude and position readers devised for and 
used with a magnifying film viewer. Time histories indicating the atti- 
tude of the model in space (angles 0 e , 0, and 0 e ), angular velocities 

(p, r, and linear velocities (u, v, w, Vp, Vy, and V), 

and angles of attack and sideslip were calculated and plotted by using the 
space coordinates and attitude readings as well as tunnel -velocity 
readings obtained from an airspeed recorder mounted in view of one 
camera. Also, by means of several spin-geometry relationships, 
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time histories were calculated and plotted for such factors as the incli- 
nations of the axis of resultant rotation and of the resultant wind from 
vertical jangles and and from each other (angle cr), and of the 

inclinations of the axis of resultant rotation with respect to the air- 
plane XZ and XY body planes jangles 5^ and A plot indicating 

the number of turns the model had made in space at any time during the 
19 seconds (full-scale time) covered by the spin and recovery was also 
made. All the time-history plots are presented in terms of full-scale 
values . 


RESULTS AND DISCUSSION 


The time-history curves of the attitude and velocity variables 
during the spin and recovery of the model are presented in figures 4 
to 13 , and motion-picture strips of the spin and recovery as filmed 
by the side and bottom cameras are presented in figures l4 and 15 , respec- 
tively. In general, the spin and recovery motion are discussed separately. 


As noted previously, the time histories were obtained for about 
3^ turns of a medium-attitude developed spin with moderate oscillations. 


(See figs. 13 to 15 •) The average values of 0 e and 0 portraying the 

attitude of the model with respect to horizontal (fig. 4) were about --43° 
(nose below horizontal) and -3° (left wing below horizontal in spin to 
pilot ! s right), respectively. The amplitude of the oscillations in 0 e 

was about ±7° and that in 0 about ±7°. The oscillatory motion of the 


spin was not completely regular, inasmuch as in any given cycle the period 
and amplitude were not exactly the same as in the preceding or following 
cycles. Along with the oscillations in 0 e and 0, there occurred cor- 
responding changes in the angles of attack and sideslip (fig. 5)- The 
average angle of attack was about 48°. The sideslip was generally nega- 
tive (outward sideslip, wind coming from left side of plane of symmetry 
during spin to pilot f s right) and averaged about -7° or -8°. The oscil- 
lations in 3 were in phase with those of angle 0 and of about the 
same magnitude. 


The angular velocity of the model about a vertical reference axis \|r e 

during the developed spin was about 90° per second (fig. 7)» This value 
increased slightly a few seconds before rudder reversal. The increase 
apparently was due to changes in the tunnel airstream associated with 
the usual practice of lowering the model a few feet just before attempting 
recovery in order to reduce materially the chances of the model crashing 
into the tunnel wall above the lower safety net. The oscillations noted 
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earlier in the attitude angles of the spin are reflected in the curves 
for p, q, and r (fig. 8) and in the remaining time-history plots. 

For example, as may be seen in figure 9> the angle between the axis of 
resultant rotation and the relative wind varied from almost zero to 
about 20°. The angle between the axis of resultant rotation and vertical 
varied from 0° to 17° (fig* 10). The relative wind direction remained 
within about 10° of vertical (fig. 10); thus, the model center of gravity 
moved almost directly toward the ground during the spin. The practically 
earthward direction of the path of the center of gravity of the model 
during the spin is also shown by the similarity of the curves for Vy 

and Vpj (fig. 11) up to the time of rudder reversal. 

As noted previously, a l~tum recovery from the spin was obtained 

by rudder reversal. (See figs. 13 to 15- ) The recovery motion required 
about 5 seconds (full scale). At the beginning of the recovery motion, 
soon after rudder reversal, an increase became evident in the oscillations 
of the attitude variables of the spin (figs. 4 and 5). By the time 
approximately 2 seconds had elapsed (l4.7 seconds on time scale) in the 
recovery attempt, the left wing had moved downward from horizontal far 
enough so that the outward sideslip reached 30°; this sideslip is 10° 
more outward than the sideslip reached during the developed spin. At the 
same time the angle of attack was about 55° • During the next second or 
so (see fig. 4, 15.8 seconds on time scale), the left wing rolled upward 
to about 20° above horizontal and the model reached a high inward side- 
slip of 16° (fig. 5)« At the time this upward movement of the left wing 
was occurring, the nose was pitching downward in space, and the angle of 
attack was decreasing to about 38°. The attainment of the 16° sideslip 
represented the first time during the whole spin and recovery that any 
appreciable inward sideslip was present. Just before the occurrence of 
the large inward sideslip, the mean rolling velocity p of the model 
doubled as the oscillations built up; just after the occurrence of the 
high inward sideslip, a momentary decrease to almost zero occurred in p. 
With the momentary decrease in p, the nose of the model began to . 
rotate downward (figs. 4 and 5)> and the angular velocities r and i|r e 

began to decrease consistently. When the model reached the stalling 
angle (a = l6°, based on unpublished results of static-force tests on a 
model at low Reynolds number) about 2 seconds after the high inward 
sideslip was obtained, the values of r and had decreased appre- 

ciably and during the next ^ second or so, they reached zero. The 

recovery from the spin appeared to have been completed at a time of about 
5 seconds after rudder reversal when the angle of attack decreased to 
below the stall. From examination of figure 11, it appears that no 
appreciable changes occurred in the resultant velocity Vp and rate of 

descent Vy during most of the recovery beyond those changes which were 
also present in the developed spin of the model in the tunnel and which 
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represent relatively small changes in tunnel velocity necessary to support 
the model. Near the end of the recovery Vy decreased somewhat whereas 

V R did not; this result would be expected since the model appeared to be 
pulling out into a glide (figs. l4 and 15 ) as the recovery was completed. 

As an aid in further discussions of the results obtained, it is con- 
venient to consider the spin and recovery in terms of the equations of 
motion about and along the three body axes. These equations are as 
follows : 


. Vp 2 pSt 

■ ~r~ c i * 0? ' I z) qr + rI xz + pqI xz 
2 

• Vr P^b / o o\ 

W* = — C m + (I z - ^ r P - (p - r 

. Vp^pSb 

I z r = — C n + (I x - I Y )pq + pI X Z " d rI XZ 


( 1 ) 

( 2 ) 

( 3 ) 


* 

U ~ 2pb 

c x - 

■ g 

sin 

e e + 

vr ■ 

- wq 

; - Zsf 




CD 

fl> 


j 


V 2|ib 

c y 

+ g 

COS 

sin f 

Oq + 

wp - wr 

„ 2 








w = VR 

2[ib 

c z 

+ g 

COS 

6 e 

cos f 

K + 

uq - vp 


W 

( 5 ) 

( 6 ) 


When p decreased to almost zero just after the large inward side- 
slip occurred, a corresponding decrease occurred in the gyrodynamic 
pitching moment (second term on right-hand side of eq. (2)); this 
decrease apparently had a large effect on the termination of the 
spinning motion. This effect appears logical because the gyrodynamic 
pitching moment is the only moment which continuously tends to hold a 
model or an airplane in a nose-up spinning attitude in balance to the 
always negative aerodynamic pitching moment* The rolling and yawing 
gyrodynamic moments are alternately positive and negative during spins 
and change with q as it oscillates above and below zero (see fig. 8 
and eqs. (l) and (3)), and the aerodynamic rolling and yawing moments are 
likewise alternately positive and negative as they seek to balance the 
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gyrodynamic moments. The products of inertia (eqs. (l) to (3)) are 
relatively insignificant since body axes and principal axes almost 
coincide . 

The large decrease which occurred in p just after the large inward 
sideslip was attained could have been due primarily to a relatively large 
negative aerodynamic rolling moment which was probably present. That such 
a rolling moment was present has been indicated by an examination of spe- 
cific numerical values of the various terms in equation (l) at the time 
the high inward sideslip was present. Making this examination required 
graphical differentiation of the time-history curve for p to obtain a 
value for the acceleration term. The value of the acceleration term was 
of limited accuracy because this differentiation was the second of two 
graphical differentiations necessary to obtain acceleration data from 
angular displacement data; however, the results can be taken as a quali- 
tative indication that the magnitude of the negative aerodynamic rolling 
moment when the high inward sideslip was present was appreciably large. 

The gyrodynamic rolling moment present at the same time was also negative 
(q, positive; fig. 9 } time, 15*7 sec). This negative gyrodynamic rolling 
moment was apparently a secondary factor also contributing to the large 
decrease in p and hence to the large decrease in the gyrodynamic 
pitching moment. 

The probable existence of a large negative aerodynamic rolling 
moment for the subject model at the time (15*7 seconds on time scale) 
a was 38°; P was l6°, and p was large is indicated by unpublished 
results of static rolling-moment measurements and damping -in- roll tests 
of various sweptback-wing configurations similar to the subject model. 

The static and damping tests indicated stabilizing dihedral effect and 
roll damping. For other wing configurations, such as highly sweptback, 
delta, or unswept wings, the values of dihedral effectiveness and of 
roll damping may be very different from those for the subject configura- 
tion. Therefore, the recovery motion for other configurations may not 
occur in a manner similar to that of the subject model. It is of inter- 
est, however, to note that, for one unswept -wing model with mass loading 
generally similar to that of the subject model, the apparent importance 
of a large amount of inward sideslip on a recovery obtained by rudder 
reversal has been noted in a recent investigation (ref. 6) in which a 
recovery motion was calculated by using rotary -balance data. 


CONCLUDING REMARKS 


A technique for determining time histories of the attitudes and 
velocities of free -spinning-tunnel models during spins and recoveries 
has been devised and presented. The technique utilizes film taken by 
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two simultaneously operating motion-pi cture cameras, one mounted at the 
side and one at the bottom of the vertical tunnel. As an initial appli- 
cation of the technique, time histories have been determined and pre- 
sented for one medium-attitude, moderately oscillatory developed spin 
and the recovery therefrom for a model representative of a contemporary 
fighter airplane. 

The results of the investigation indicated that the oscillatory 
motion of the spin was not completely regular inasmuch as the period 
and amplitude of any one cycle of the time-history curves were not 
exactly the same as those for the preceding and following cycles. 

After rudder reversal, increases occurred in the oscillations, partic- 
ularly in the angles of wing tilt and sideslip and the recovery was 

completed about 1^ turns after the rudder was reversed. One factor which 

apparently had a significant part in the recovery motion obtained in this 
investigation was a relatively high inward sideslip reached during the 
oscillations. This sideslip apparently was the primary cause of a sudden 
large decrease in rolling velocity; this decrease caused a corresponding 
large decrease in the gyrodynamic pitching moment which had been holding 
the model in a nose-up spinning attitude. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., November 7, 1955* 
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APPENDIX 

TECHNIQUE OF DETERMINING TIME HISTORIES OF ATTITUDES AND 
VELOCITIES OF MODEL DURING SPIN AND RECOVERY MOTIONS 


The motion-picture camera at the bottom of the tunnel was a 35 - 
millimeter camera with a wide-angle lens. This camera was aimed verti- 
cally upward along the center line of the tunnel and was not moved from 
this attitude as it photographed the motion of the model. The camera 
at the side of the tunnel was a l6-millimeter camera which was yawed 
and pitched as necessary to follow the motion of the model. The two 
cameras were not mechanically synchronized but both recorded accurately 
the passage of time by means of timing devices located in their views. 

From film from the side camera, readings of the angle rj for four 
points on the model for numerous times were obtained. These four points 
are the nose and tail ends . of the X body axis and the right and left 
wing tips. The angle t} readings were obtained by making use of lines 
painted on the tunnel -wall panels behind the model. (See sketch, fig. 16.) 
In making the trigonometric calculations necessary for laying out points 
through which these r| lines (fig. l6) were painted on the walls, con- 
sideration was given to the dimensions of the tunnel (dodecagonal cross 
section) and to the distances of the side-camera lens from the reference 
horizontal plane P^ and from a vertical reference line through point Q 

when the camera was pitched t|° and yawed arbitrarily. These latter 
distances varied with t) in accordance with the relationships 

Sjj c = 0.1544 cos n - 0.7500 sin t| ( 7 ) 

and 

Sy = 0.1544 sin t) - 0.7500 (l - cos T]) (8) 

V ^ c 


(See sketch, fig. IT*) 

From film from the bottom camera, readings at numerous instants were 
obtained for the angles A and e for the same aforementioned four 
points on the model for which r\ was read. Readings were also obtained 
from the film from the bottom camera for the angles \|r e , v, and £ and 

the distance p «. The sketch in figure l8 illustrates these angles and 

n, 1 

distances as they appear on the film. Special clear plastic templates 
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were designed and made to enable the reading of these factors on the film 
with an available magnifying film viewer. 

The Tj readings from the film from the side camera were synchronized 
with the film readings from the bottom camera timewise by plotting the tj 
readings against time and picking off values for the same instants at 
which the bottom-camera film readings were made. Then, values of Sy 

for the four aforementioned points on the model were calculated by use of 
the equation 


s 


V 


( 10.95 - 23 .48 tan A - s H c cos e)tan tj + s v^c cos £ 
cos e + tan A tan tj 


( 9 ) 


which is derived as follows: Reference to the sketch in figure 19 

indicates that 


tan tj 


S _V - s V,c 
S H ' s H,c 


and 


tan A 


10.95 - s R f 

23 .48 + Sy 


( 10 ) 


(H) 


From equation (10), 


S V " s V,c . 

Stt = T + s„ 

H tan r) c 


and, as may be seen from figure 19 , 


3 H 


' = Sg cos e = 


tan tj H 




cos e 


From equation (ll). 


( 12 ) 


Sjj* = 10.95 - (23.^ + Sy^tan A (13) 

Equating the right-hand sides of equations (12) and (13), expanding, 
and rearranging gives equation ( 9 ). 
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The attitude angles of the model in space were calculated by the 
relationships 


e e = 


sm 


-1 Vjt ~ s V»n 


(14) 


0 


sin 


-1 


s V,w l ~ s V,wr 
h 


(15) 


and 


0e 


sm 


-1 


sin 0 
cos 0 e 


(16) 


(See fig. 20.) ^It is of interest to note here that, if a model were 
to he in an inverted attitude rather than in an erect attitude, equa- 


tion (15) would take the form 0 = sin“l wr - ^ The angular 
velocities p, q, and r were calculated from 


P = 0 e " sin 6 e 


• • 

q = 0 e cos 0 e + t e cos 0 e sin 0 € 


(17) 

(18) 


and 


. . 

r = -0 e sin 0 e + \|r e cos 0 e cos 0 e (19) 

(angular velocities about body axes in terms of Eulerean-type space angles 

• . . 

similar to those of ref. 7). Values of 0 e , 0 e , and i|r e were obtained 

graphically from 0 e , 0 e , and \|/ e . Values of the resultant angular 

velocity Cl were obtained by vectorial addition of the mutually per- 
pendicular angular velocities p, q, and r. 

The linear velocities u, v, and w were calculated by the fol- 
lowing equations: 


dp 


u = ± 


H _ 

— cos v cos 0 P + -fp^. sm v 
dt e dt H 


cos 0 C 


V. r sin 0 P 
V 


( 20 ) 
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dp t u 

v = ^2 sin(v ± B)cos 0 ± p^. cos(v ± B)cos 0 + Vy sin 0 


( 21 ) 


cLp_j n 6» 

w = -f — — cos(v + A ± B)cos ob - p sin(v + A ± B)cos a + 
at 3 at H j 


V v cos 9 e cos 0, 


( 22 ) 


When a choice of signs is indicated in equations (20), (21), and (22), 
the upper sign applies if v = and the lower sign applies if v = v n . 

As an aid in visualizing the angles used in equations (20), (21), and (22) 
and the significance of the relationships, refer to the sketches in fig- 
ure 21. The instructions regarding signs would be different from those 
used if 0 e were to become positive or if the model were to be in an 

inverted attitude. 


Values of for use in equations (20), (21), and (22) were 

obtained by multiplying the tangent of the angle represented through 
distance p on film-magnifier glass (fig. l8) times the factor 

^23-48 - Sy The latter factor represents the vertical distance of 


the center of gravity of the model above the lens of the bottom camera. 
Values of s„ were determined simply from knowing s. r , s , , 

V y eg V,n V y "0 

the fuselage length, and the location of the model center of gravity 

3-Ptj df 

along the fuselage length. Values of and — were obtained 

dt dt 

graphically from and from £. Values of a A, and B were 

calculated from the relationships (see figs. 20 and 22) 


sin = cos 9 e cos 0 e 

sin A = tan 0 tan 
cos 0 e 


cos B = 


cos 0 


Values of Vy were obtained from 


ds T 


V = y + Y 

v V v tunnel ^ 


(23) 

(24) 

(25) 


(26) 


where the latter term was determined graphically. 
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Values of the resultant linear velocity Vp were obtained by vec- 
torial addition of the mutually perpendicular linear velocities * u, v, 
and w. Angles of attack and sideslip were calculated from 


a = tan“^ — 
u 


(27) 


and 

3 = sin" 1 X (28) 

K 

The angles between the axis of resultant rotation and the XZ and 
XY body planes were calculated by 


and 




'ft 


sin 


-1 


r 

ft 


(29) 


(30) 


The angle between the axis of resultant rotation and the direction of 
the resultant wind was calculated by 

cos a = cos cos Typ cos^5^ - &vr) + si - n sl - n t vr (31) 

(See fig. 23.) The velocity along the axis of resultant rotation was 
calculated by 


V = cos a 


(32) 


The inclination of the axis of resultant rotation from the vertical was 
calculated by 


cos K 


ft 


cos oc^' 


cos a^ + 


sin 


“1 


sm 


0 ^ 


cos cos 


-1 


± cos 


sin ay 


P + q 2 


sin 




(33) 


which is derived fully in reference 6. Equation (33) represents a spher- 
ical trigonometric method for determining and makes use of known 

values for the attitude of the airplane with respect to the vertical and the 
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inclination of the axis of resultant rotation with respect to the 
airplane. In the present computation, the negative sign in the last 
term of equation (33) was used when 5^ and 0 were of like sign, and 

the positive sign was used when 6^ and 0 were of opposite sign. The 

inclination of the resultant linear velocity vector from vertical was 
calculated from 


COS KyR = -X (34) 

V R 

The number of turns the model had made in space at any time during the 
spin and recovery was determined directly from the values of \| / read 

from the film. 

It should he noted that extreme care was used in establishing the 
accuracy with which the film could he read and with which slopes of 
curves could be taken by graphical procedures. These accuracies were 
used judiciously in fairing the measured data and the subsequent calcu- 
lated results. 
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TABLE I 


DIMENSIONAL CHARACTERISTICS OF AIRPLANE REPRESENTED BY MODEL 


Overall length, ft 


37-5 


Wing: 

Span, ft 

Area, sq ft 

Sweep! a ck at quarter chord, deg 

Dihedral, deg 

Incidence at root, deg 

Incidence at tip, deg 

Section at root 

Section at tip 

Aspect ratio 

Mean aerodynamic chord, in 

Leading edge of c rearward of leading edge of 
root chord, in 


.... 37.1 

.... 287.9 

.... 35 

.... 3 

.... 1 

.... -1 

NACA 0012-64 modified 
NACA 0011-64 modified 

.... 4.8 

.... 97.0 

46.0 


Ailerons : 

Area, sq ft 39*2 

Span, percent semispan 50.1 

Hinge-line location, percent chord 70 


Horizontal tail: 

Total area, sq ft 35-0 

Span, ft 12.8 

Sweepback at quarter chord, deg 35 

Elevator area rearward of hinge line, sq ft 10.1 

Distance from center of gravity to elevator hinge 

line at root, ft l8.0 

Dihedral, deg 10 

Incidence, deg 0 

Vertical tail: 

Total area, sq ft 33-7 

Sweepback at quarter chord, deg 35 

Rudder area rearward of hinge line, sq ft 8.7 

Distance from center of gravity to rudder hinge 

line at root, ft 17*3 

Unshielded rudder volume coefficient 0.0143 

Tail-damping ratio 0.0204 

Tail-damping power factor 0.000290 


Ol|N O l| X 
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TABLE II.- MASS CHARACTERISTICS FOR MODEL 

[Values are given in terms of full-scale airplane; moments of 
inertia are about center of gravity] 


Weight, lb 

p at an altitude of 15,000 ft 


1^, slug -ft 2 

l y, slug-ft^ 

l z , slug -ft 2 

^ ^ 

mb 2 



mb 2 


13,103 

25.^5 

0.221 

0.133 

9,255 

25,859 

33,205 

-296 X 10“ 4 
-131 x 10" 4 
427 x 10 -4 
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(a) 0 e and \|r e = 0. 




Z 

(c) 9 e and \(r e = 0, and in this case 0 = 0 e . 
Figure 1.- Body system of axes and related angles. 
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Figure 2.- Three -view drawing of the model. 
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Figure 3*- Photograph of the model. 



Figure 4.- Variation of attitude angles 0, 0 e , and 0 e with time during 

spin and recovery of model. 


MCA TN 3611 



Figure 5*- Variation of angles of attack and sideslip with time during 

spin and recovery of model. 
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Figure 6.- Variation of linear velocities u, v, w, and Vr with time 

during spin and recovery of model. 
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Figure 7»- Variation of rate of rotation about vertical axis with time 
during spin and recovery of model. 
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Figure 8.- Variation of angular velocities p, q, r, and 0 with time 
during spin and recovery of model. 
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Figure 9-- Variation of angle between relative wind and axis of resultant 
rotation with time during spin and recovery of model. 
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Figure 10.- Variation of angles between axis of resultant rotation and 
vertical and between relative wind and vertical with time during spin 
and recovery of model. 
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Figure 11.- Variation of linear velocities V, V y , and V R with time 
during spin and recovery of model. 
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Figure 12.- Variation of angles 5^ and (locating axis of resultant 

rotation with respect to XZ and XY body planes, respectively) with 
time during spin and recovery of model. 
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Figure 13- - Variation of turns completed about vertical axis with time 
during spin and recovery of model. 
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Figure 14.- 


L-90537 

Motion-picture strip of model spin and recovery as filmed 
from side camera. 
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Frames 
26$ to 236 


Frames 

287 to 290 


mm 


m 

! 4 


Frame 

287 

model 

in 

dive 


Frame 271 j 
about 1^ 

turns after 
reversal 
(model re- 
covers ) 


- 

'® 4 

HR* « 

3K 


Figure 14.- Concluded. 
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Frames Frames 

1 to 6 7 to 12 



Frames 
13 to 18 



frames 
19 to 2 b 



Frames 
25 to 30 



L-90533 

Figure 15 •- Motion-picture strip of model spin and recovery as filmed 

from bottom camera. 
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Frames Frames 

31 to 36 37 to 42 



Frames Frames 

43 to 48 49 to 54 



Frames 

55 to 60 



Figure 15*- Continued. 
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Frames 
6l to 66 



Frames Frames 

67 to 72 73 to 78 



Frames 
79 to 84 



Figure 15 . - Continued. 


Frames 

85 to 90 
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Frames 

91 to 96 



Frames 
97 to 102 


Frames 

103 to 108 


Frames 
109 to 114 



Figure 15. - Continued. 


Frames 
115 to 120 
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Frames 
127 to 132 



Frames Frames 

133 to 138 139 to 144 



Figure 15. - Continued. 


Frames 
145 to 150 
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Frames 
151 to 156 


Frames 

157 to 162 


Frames 

163 to 168' 



I Frame 139} 

I about li turns 

latter rudder 
[reversal 
(model recovers)! 



Frame 

model 

dive 


Figure 15.- Concluded. 


Frames 

169 to 171 



L-91698 
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Figure 16 .- Sketch indicating model and lines on tunnel walls as seen 

from film taken by side camera. 
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Figure 18 .- Sketch illustrating model, side camera, and pertinent angles 
as seen from film taken by bottom camera. 
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Figure 19- - Sketch indicating tunnel, reference planes, angles, and dis- 
tances used in determining Sy of a point on the model. 


Horizontal reference 
plane P w 


Center of bottom camera lens 


Reference point Q 


Point on 
model 




23.48 


' - Center of side 
camera lens 
( arbitrary camera 
position ) 


Vertical reference 
plane Py 



50 


NACA TN 3611 



Vertical plane through Y body axis 
Vertical plane through Z body axis 

I | YZ body plane 
| | Vertical plane through X body axis 

Figure 20.- Sketch showing spherical trigometric relationships between 

sin (f> 


angles 0, 9 e , , and a^. Sin = 


cos 0 e 
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dA 

if 


(a) 0 positive; v n positive. 



(b) 0 positive; v-t positive. 


Figure 21.- Vector diagrams indicating various relative attitudes of pos- 
sible linear velocity components in horizontal plane. Horizontal pro- 
jections of body axes shown as viewed from below model, for cases (a), 
(k)> ( c )> a-nd. (d) . Origin of vector system is center of gravity of 


model. 


Vector direction of points toward various possible loca- 


tions of tunnel center line, as seen on a frame of film from the bottom 
camera. Applicable v selected in accordance with definition in 
symbols. In spin and recovery, 0 e is always negative. 
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(c) 0 negative; v n positive. 



(d) 0 negative; v t positive. 


Figure 21.- Concluded. 


Vertical direction 


CO 



Figure 22.- Sketch indicating spherical trigonometry relationships per 
tinent to determination of angles A and B. 
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NACA - Langley Field, Va. 



'Ol 

-f=" 


cosBT = cos(8 ^-S vr ) cos t vr 
sin t vr 


sin 3E 


sin EE 


cos <7 * cos T^cosEZ+sinr^sinlSE cos(90-3E) 
cos <7 *cosTqCost vr cos( 8^-8 vr ) +sin r, 


Sin T wn 

m sin3E 


sin IE 

cos < 7 s cosT^cosTypCos( 8 ^- 8 y R ) + sin T^sin Ty R 

Figure 23 .- Sketch indicating vector diagrams and spherical trigonometry- 
relationships pertinent to determination of angle 0 . 
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